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O ATMOSPHERIC DYNAMICS OF RED SUPERGIANT STARS 
<N AND INTERFEROMETRY 
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' ' . Abstract. We developed a 3D pure LTE radiative transfer code to de- 

' rive observables expected for RSGs, with emphasis on small scale struc- 

, tures, from radiative-hydrodynamic (RHD) simulations of red super- 

Ph- giant stars (RSGs) carried out with C05B0LD ( |Freytag et al., 2002D . 

JL ' We show that the convection-related surface structures are observable 

< I , with today's interferometers. Moreover, the RHD simulations are a 

■ great improvement over parametric models for the interpretation of 
^ I interferometric observations. 



^ 1 Introduction 
ON 

Red supergiant (RSG) stars are irregular, small- amplitude variables and they 
CO ■ have a very complex spectrum where strong molecular absorption by TiO and 
other molecules dominates. This leads to an ill-defined continuum. Furthermore, 
■ Josselin fc Plez (2QQ7)| find that high-resolution spectroscopic time-series observa- 



o 



I tions of a sample of RSG reveal that the variable velocity field in the atmosphere 
■ probably has convective origin. The velocities are supersonic and vary with a time- 
scale of a few 100 days and the resulting spectral lines have a strong broadening 
together with asymmetries and wavelength shifts. For all these reasons, spectral 
synthesis of these objets is really difficult to achieve. 

The observed strong line asymmetries indicate that convection consists of giant 
cells, as suggested by radiative hydrodynamic modeling (see Freytag contribution 
in this volume). 

1 . 1 Radiative-hydrodynamic simulations 

The numerical simulations analyzed here are performed with "C05B0LD" ("Con- 
servative code for the Computation of COmpressible COnvection in a BOx of L 
Dimensions, 1=2,3") developed by Freytag, Steffen, Ludwig and collaborators (see 
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[Freytag et al., 2002] ). The RSG simulations presented employ the global star-in-a- 
box setup: the computational domain is a cube, and the grid is equidistant in all 
directions. Radiation transport is strictly LTE. The grey Rosseland mean opacity 
is a function of pressure and temperature. The necessary values are found by 
interpolation in a 2D table which has been merged at around 12 000 K from high- 
temperature OPAL data ( Iglesias et al., 1992 ) and low-temperature PHOENIX 
data ( Hauschildt et al., 1997D by Hans-Giinter Ludwig. 

Some more technical informations can be found in the C05B0LD Online User 
Manual and in the contribution by Freytag (this volume). 



1.2 A posteriori radiative transfer in the hydrodynamical simulations 

We performed 3D pure LTE radiative transfer calculations in snapshots of the 
3D hydrodynamical simulations (the RHD model presented in this work has a 
solar composition, 235^ grid points, M = I2M0, L = lOOOOOL©, Tg// = 3482K, 
log^ = —0.39, and R = SSOi?©), taking into account the Doppler shifts caused 
by convective motions. The emerging intensity at the top of one ray of the RHD 
computational box, is computed by summing the contribution of the source func- 
tion at different depths. 

To reduce the computing time, the extinction coefficients per unit mass are pre- 
tabulated as a function of temperature, density and wavelength with a resolution 
of 0.01 A using MARCS ( |Plez et al.,T992l [Gustafsson et al., 2003^ opacity data. 



We checked that this resolution is sufficient to ensure an accurate characteriza- 
tion of line profiles even after interpolation at the Doppler shifted wavelengths 
dChiavassa et al., 2006] ). 

The a posteriori radiation transfer treatment is used to compute : (i) spectra 
at high and low (Opacity Sampling) spectral resolution that are compared with 
spectroscopic and spectrophotometry observations; (ii) intensity maps to compare 
with inter ferometric observations. 

Concerning the spectral synthesis. Fig. [T] shows that RHD models approximately 
reproduce the 12 fim H2O line width without the need for micro- or macrotubu- 
lence (which is needed in the one-dimensional model). The shifts of the calculated 
line center span values between -1 and -3 km/s with respect to the hydrostatic 
model. The profiles are asymmetric. 



2 Interferometric imaging of red supergiant stars 

The characterization of stars beyond simple limb-darkened circular disks is a very 
important task to be achieved. Our 3D radiative transfer code is employed to 
derive observables expected for RSGs with emphasis to small scale structures. 
We produced a time- sequence, with a time-step of about 23 days covering about 
one stellar year, of maps of monochromatic intensity at two different wavelengths. 
These wavelengths probe different depths in the stellar atmosphere since they 
correspond to the continuum at 1.6/im, where the H~ minimum opacity occurs 
and one can see deeper in the photosphere, and a nearby CO line. Two examples 
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Fig. 1. The time sequence computed for the RHD model described in the text (time-step 
of about 23 days) is compared to observations (dots; Ryde et al., 2006) the sharp decrease 
below 818.38cm~^ is an instrumental artifact). The dashed line is a one-dimensional 
hydrostatic MARCS model calculation convolved with an exponential profile of 13km/s 
width. 



of the maps are displayed in Fig. [21 together with the radially averaged intensity 
distribution. Absolute model dimensions have been scaled to approximately match 
the angular diameter of 35mas. The contrast between large bright granules visible 
in the continuum (upper- left panel), and the smaller structures seen in the CO 
line (bottom-left panel) is striking. 

We computed the visibilities curves from a large amount of intensity maps for 
different snapshots and rotated images (this analysis is extendable to other mod- 
els and wavelengths). The visibilities carry a clear signature of deviations from 
circular symmetry (Fig. [3]) because the small scale structures lead to a significant 
scatter. The emerging intensity at a given wavelength at a given position on the 
star depends on the opacity run through the atmosphere below that point. This 
in turn depends on the temperature and density stratification. Observations at 
a wavelength in a spectral line and in the nearby continuum probe thus different 
depth and thus different temperatures in the photosphere. In order to characterize 
the convective pattern, high spectral resolution is needed to determine the varia- 
tions of the interferometric data between different spectral features, and between 
features and continuum. 

Furthermore, the visibility strength at 100m is 0.02 — 0.03, much higher than what 
we expect for a uniform disk distribution. This is due to the small wiggle struc- 
tures visible in the intensity maps (Fig. [2]) that are not spatially resolved. 
The convective-related surface structures are then observable with today's inter- 
ferometers. 
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Fig. 2. Logarithmic maps of monochromatic intensity (a grey scale is used and lighter 
shades correspond to higher intensities) and scatter-plot of the radially averaged intensity 
computed for 48 snapshots from RHD simulations described in the text with a time-step 
of about 23 days (grey; the black line is the averaged profile). Upper row : continuum at 
1.6/im (H~ minimum opacity). Bottom row : Nearby CO line. 

2.1 Radiative-hydrodynamic models and Interferometry 

RHD simulations are necessary for a proper, qualitative and quantitative analysis 
of the surface of cool stars, in order to have parameter-free estimates of convection. 
RHD models are ab initio^ time-dependant, multi-dimensional, non-local and they 
contain as much physics as possible. The actual limitation is the CPU-time. RHD 
models are a great improvement over parametric models for the interpretation 
of the observations but they have to be validated with other techniques such as 
spectrometry and spectrophotometry ( [Chiava ssa et a l., 2006 ). 
In particular, VSfl will produce reference observations for RSG stars from where 
the typical size of surface inhomogeneities, their contrast (in temperature and flux) 
and the order of magnitude for the time-scales of variations can be derived. These 
quantities are relevant for validating present models and for pushing further model 



-•^VLTi Spectro-Imager is proposed as second generation VLTI instrument providing the ESO 
community with the capabihty of performing image synthesis at miUi-arcsecond angular resolu- 
tion 




Fig. 3. Scatter-plot of visibilities (grey) as a function of projected baseline length ori- 
entation, and the average profile (black), for 48 snapshots and rotated images computed 
from the RHD simulation described in the text. The dashed curve is the visibility of 
a uniform disk with a diameter of 35mas. Left column : continuum at 1.6/jm. Right 
column : Nearby CO line. Bottom panels are enlargements of the upper ones. 



developments. 

2.2 Future developments of the radiation transfer code 

The main characteristics of the 3D radiative transfer code we developed have 
been reported in Sect. II. 2[ On a short time-scale, future developments consist in 
introducing scattering and in achieving a 3D polarized radiative transfer code. This 
is particularly important to properly characterize the convective motions of RSG 
by constraining the degree of asymmetries and associated temperature contrasts. 
For this purpose, we plan spectropolarimetric observations. On a long time-scale, 
Non-LTE radiative transfer is required in order to properly analyse RSG spectra, 
since our final objective includes reliable determination of abundances for these 
stars. 



3 Conclusion 



The comparison of our RSG RHD models with interferometric and spectroscopic 
observations is now possible with our radiative transfer code. 
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We found that convection-related surface structures are indeed observable and the 
interferometric visibilities carry an evident signature of deviation from circular 
symmetry. The phase and related closure phase show clear deviations from the 
axisymmetric case as well. The statistical approach applied in this work will be 
used to set a range for visibility /phase variations due to the photospheric convec- 
tion to predict further observations (Chiavassa et al., in preparation). 
There are few available interferometric observations suitable to constrain the mod- 
els in terms of visibility and phase. Nevertheless, since RSG are prime targets for 
Interferometry, thanks to their large diameter and to their high-peak infrared lu- 
minosity, a lot of observations will be available in the near future. In this context, 
we think that a theoretical and observational simultaneous approach is really im- 
portant to understand the atmospheric dynamics. 

We thank Bertrand Plez and Eric Josselin for the useful and enhghtening discussions and Bernd 
Freytag for the RHD models. We thank the "Centre Informatique National de I'Enseignement 
Superieur" (CINES) for providing us with computational resources required for part of this work. 
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